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Abstract 

We present a dark energy model with a double exponential potential 
for a minimally coupled scalar field, which allows general exact integration 
of the cosmological equations. The solution can perfectly emulate a A- 
term model at low-medium redshift, exhibits tracking behavior, and does 
not show eternal acceleration. 

1 Introduction 

Among the many problems posed by the discovery of positive acceleration of 
the Universe, one of the most puzzling is that of eternal acceleration. In the 
case of a true cosmological constant, as well as for most of the models proposed 
in the literature, the present time accelerating era seems doomed to last forever. 
This creates problems in string theory, where asymptotic flatness of spacetime 
is required 0. Several possible solutions are based on cycles in the evolution 
of the Universe 0, or on the introduction of negative potentials for the scalar 
field 3 . A much simpler approach is obtained with a negative cosmological 
constant, together with a suitable potential 0]. However, this last approach 
poses the same type of problems which make unsatisfactory the simple use of 
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a standard positive A. A common feature of many of these models is that 
they foresee a recollapsing universe in the context of a spatially flat metric. 
But a recollapsing universe can also be obtained simply by setting k = +1, 
an option which has not been ruled out by the BOOMERANG and MAXIMA 
missions, as instead people often say. On the contrary, the most recent data 
by Wilkinson Microwave Anisotropy Probe (WMAP) have reopened the way to 
nonflat cosmological models 0, [H]. 

In this paper, we remain in the context of a flat space (k — 0) and propose 
a very simple and attractive model with a minimally coupled scalar field plus 
dust. We shall find a general exact solution of the equations, and show that it 
meets all the requests which are usually made for a dark energy model, including 
tracking behavior and switch off of the acceleration. Moreover, we shall not get 
a recollapsing universe. 

Let us consider the following potential (in units such that 8irG = 1): 



This potential has already been discussed in the literature (see Refs. [7], [8], 
for instance). In Ref. which is a paper closely related to this one, the 
potential is not equal but only similar to that in Eq. Comments on the 

attractor properties of the (generic) solutions due to such a kind of potential 
in connection with initial conditions and the values of model parameters have 
also been given, without trying to find any exact solution [7|. On the other 
hand, Ref. P inverts our point of view, since it postulates a special scale factor 
behavior, so deducing a potential form similar to the one in Eq. Qvith some 
relevant differences with respect to it. (The main point is the different relation 
among the values of the potential parameters.) The potential also turns out 
to be a particular case of the wider class discussed in Ref. 0] . For a discussion 
on the possibility of finding a general exact solution for the kind of potentials 
to which the one in Eq. (P) belongs, however, see Ref. [TU| . 

In any case, here we do not present a new technique for finding solutions; 
instead, we select the particular class of potentials in Eq. JQl and a treatment 
of the related solutions, which allows us to obtain the desired results and puts 
them in greatest evidence. 

Some comments are now needed on Eq. 

(1) The second form of the potential can also be interpreted as containing a 
negative cosmological constant, fine tuned with the parameters of the potential 
itself. The first form shows that this is in fact only a matter of interpretation. 
Such a form as a whole, of course, still retains the right to be used as a potential 
in itself. 

(2) In particular, if we set A = B, the expansion of V(ip) around its minimum 
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V(<p) ~ -{3B 2 )p 2 + — (3B 2 )p A + (2) 

showing that it can be interpreted as a massive self-interacting potential, which 
is quite reasonable in the situation we are studying. We shall see that posing 
A = only affects the asymptotic value of tp. As a matter of fact, we consider 
A ^ B just for sake of generality. 

(3) The real problems with potential Q are the particular value of the 
exponent, and the link between the mass term and the self-interaction term 
in Eq. J3J). We are not really able to explain them in terms of fundamental 
theory 1 . It is however intriguing that, as we shall see, the dynamics entailed by 
this potential on the minisuperspace {a, p} admits a Noether symmetry. 

After presenting the cosmological solution in Sec. II, we analyze it in Sec. 
III. Finally, Sec. IV is dedicated to some conclusive remarks. 



2 Solution 



Let us consider the Einstein equations, together with the Klein-Gordon equation 
for the scalar field, in the flat case 



.a i a\ 1.2 



3-ff = P v +Pm, (3) 

2- + ( - ) +■ ^P 2 - V(<p) = 0, (4) 
a \a J 2 

<p + 3Htp + V'(<p) = 0. (5) 

Let us concentrate on the last two equations and put aside the first one, 
which is, as is well known, a first integral. The reason for such a choice is that 
these equations can be derived by the Lagrangian 

L = 3aa 2 - ^a 3 ip 2 + a 3 V(<p). (6) 

In close analogy with Ref. we now set a transformation of variables, 
which is however not the same as there 



a 



■ 3(u 2 -v 2 ) 2, u + v 

3 = A 2 , P=\ ^Og 

4w V o u — v 



^=\/o lo g \-Poo, (7) 



where uj 2 = SAB and p^ = y/ 2 / Z\og(B / A) . With this change, we get for L 

L = \{ii 2 -v 2 ) +v 2 , (8) 



But, about this, see Refs. EU,C2| for some dimensional considerations on exponential 
potential, once one wants power law solutions for the scale factor. 
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and we see that the variable u is cyclic, so that we have got a Noether conserved 
quantity 

dL 2 . 3 



._> ■ , 2 y/"V(<p)(a + a<p). (9) 



L0 Z LO* 



We can now easily obtain a general exact solution 

u = u\t + «2 > v = V\ sm{u>t + V2). (10) 

Thus, we have four integration constants (u\ t U2, V\, 1)2) and two free pa- 
rameters (lo and ifoo , which are in turn derived from A and B in the potential) . 

First of all, let us set u 2 = v 2 = 0. This is more than needed to fix a time 
origin, i.e., a(0) = 0, but it is important to select the tracking solution. At the 
end of Sec. Ill we shall consider a weaker condition. 

The second step consists in fixing the time scale; we pose, then, to = 1. 

Finally, let us denote by Ho the present value of the Hubble parameter [but 
measured in units of (age of the universe) -1 ], which in our units is of order 1, 
even if it does not coincide with the value of the usual h. 

When all these conditions are imposed we obtain the following expressions 

3 _ (3H a - 2) esc 2 lo sin 2 Lot + (2lo cot Lo-3H )t 2 (n] 
a ~ 2(lo cot lo- 1) ' ^ ' 

up = <J~^\og j (tsinuiy/ (3Ho — 2uo cot lo)(lo cot lo — 1) 

— yj(2 — 37<o)(l — w cot w) sin Lotj j t sinLoy / (3Ho — 2lo cot lo) (lo cot lo — 1) 

+V (2 - 3H ) (1 - w cot to) sinurfj + ip^ , (12) 

2t(3Ho — 2lo cot lo) — (3Ho — 2)u> esc 2 lo sin 2ut 
3t 2 (3Ho — 2lo cot lo) — 3 (3Ho — 2) esc 2 lo sin 2 cut ' 

y, t j = 4(3H - 2)u; csc 2 cj sin 2 Lot 

3t 2 (2LO cot lo — 3Ho) + 3(3Ho — 2) esc 2 lo sin 2 Lot ' 
. 2 8 (3Ho — 2) sin lo (3Ho sin lo — 2 cos lo) (sin Lot — Lot cos ut) 2 

V 3(t 2 (LO sin 2lo - 3H sin 2 lo) + (3H - 2) sin 2 Lot) 2 ' 

From Eq. (jl 1|> it is clear that a, a, and a are all oscillating functions. 
Actually, a still oscillates to infinity, when it behaves like 

2§ (-2 + 3 H ) lo 2 csc(lo) 2 . . 

aoc i J ^ T cos(2 Lot). (16) 

3 (-1 + co cot(w)) 3 (i 2 (-3 Ho + 2 w cot (a;))) 3 

It is interesting to note that the right-hand side in Eq. Ijltifl is a damped oscil- 
lating function, unless Ho = § ; however, since Ho is not given by the usual h 
parameter here, as noted before, such a value is forbidden by the recent obser- 
vational data. 

From the equations above it is straightforward to derive the expressions for 
tt v = pyl (3H 2 ) and w v = p v / p v (p v being the pressure connected to the scalar 
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field), which are not given explicitly here, since they are really messy. It is also 
clear that (^oo only enters into ip(t), and determines its asymptotic value. We 
are, thus, left with uj and TLq as free parameters. 

3 Analysis of the solution 

First of all, let us observe that Eq. I|13|) . no matter which are the chosen values, 
shows that H(t — > oo) rs 2/ (M). That is, the behavior in the future approaches 
pure dust, and the acceleration, if ever present, is switched off. This, of course, 
does not mean very much, until we do prove that the model can match the 
observations. 

As we are not interested in best fits with real observational data, but rather 
in comparing the behavior of our model with a pure A-term model, it is sufficient 
to choose reasonable values for ui and Hq. Since Hq ~ 1, let us set this value for 
the moment and look at the plot of Q 0ip versus uj (see Fig. 1). We discover that 
large ui values yield nonphysical results. The best choice seems to be w < 2. 
Provided this limitation is set, the value of Qq v is actually not very sensible to 
uj, so that this parameter is in fact almost degenerate. 

Finally, we choose (in our units) 

H a = 0.95, uj = 0.3 =*> fl 0v > ~ 0.72, (17) 

which seems quite reasonable. Let us also consider a cosmological constant 
model with 

H = l, Q A = Q <p ra 0.76. (18) 

We have to insist on the fact that these values are just indicative. What we 
want to say is that, after considering one of the two models as fiducial, they 
are perfectly equivalent if one takes rather small differences in the free parame- 
ters. Indeed, comparing the two luminosity distance curves we obtain the strict 
overlapping in Fig. 2. 

It should be clear from this that the two models will stay undistinguishable 
not only for the present observations with SNIa, but also with near future 
improvements, like those supposed to come from Supernova/Acceleration Probe 
(SNAP) mission |15| . Only a tremendous advance in the precision could make 
the difference. 

A comparison of the two models with real data on peculiar velocities of 
the galaxies is also made in Ref. and we again find good agreement and 

impossibility to distinguish. 

A study of the consequences on the cosmic microwave background radiation 
(CMBR) spectrum is now in preparation, but a first check seems to give the 
same type of answer |T7) . 

In order to complete the discussion, let us look at the plot of the acceleration 
d(t) (Fig. 3). We see that it evolves towards alternative periods of accelerations 
and decelerations, which justifies the title of this paper. With our choice for 
the parameters, all this happens in the far future, and we now live in the first 
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period of acceleration. The reason why it is so is connected to the coincidence 
problem. In any case, we have a hint on a possible way to face the problem. 
It may be that a more sophisticated model gives a less steep damping of the 
oscillations. 

Let us now show the plots of log 10 p v versus log 10 a and w = w v versus 
lo g 10 a (see Figs. 4 and 5). After the above discussion, the first plot should be 
self-explicative: it has the shape usually obtained with qualitative arguments, 
but in this case it is obtained by exact integration, which seems noteworthy to 
us. In the second plot the present time is marked by the zero point and seems 
to be very special. We can imagine to be in the far past or in the far future, 
and make the same analysis with largely different values for Qq v . Considering 
such values in the range {0.01 ~ 0.93}, we obtain the overlap in Fig. 6, showing 
that the coincidence is only apparent. 

Last but not least, we can show that the examined solution also exhibits the 
socalled tracking behavior. 

So, let us release the condition U2 = i>2 = 0, and substitute it with the more 
general ui — Vi — x, which is the minimum necessary to have a(0) = 0. Then, 
let us also consider all the other constraints introduced in Sec. II. Allowing 
x to vary over a very narrow range around zero, Fig. 4 is modified into Fig. 
7, showing that a large set of initial conditions eventually evolves towards the 
same final behavior. Even if x is increased, the qualitative behavior does not 
change much, but the time needed for reaching the tracking solution is now 
larger, leaving still opened a sort of fine tuning problem. In any case, let us 
observe that x is only a mathematical object, without any obvious physical 
interpretation. The relevant physical quantity p v is not fine tuned, as we have 
shown. 

4 Conclusions 

The model presented here, in our opinion, shows a possible solution for all the 
main problems posed by the dark energy theory, at least from a mathematical 
point of view. 

The physical attitude should be however more careful. Though very simple 
and reasonable, the potential proposed here has at the moment not satisfactory 
justification, independent of this model. 

We have seen that, were this the real situation, it could be probably impossi- 
ble to discriminate from a A-term model with any realistic plan of observations. 2 . 
As for the switch of the acceleration in the future, it is of course unobservable 
in principle! 

The main result of this paper seems to be an indication of a possible way to 
solve the problem of switching off the acceleration by means of minimal tools. 
It also suggests that oscillations of the scalar field could be an interesting field 

2 Note that this tendency also results from the very recent WMAP data |18| 
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of study also in a context very far from inflation and reheating. (For some 
considerations about this, for example, see Refs. [T§|, [2D1 -) 

A last remark has to be addressed to the fact that the radiation dominated 
era is not recovered at early times. As a matter of fact, this is due to the 
model itself and the way it is built. Adding a radiation component is not 
simple on an analytical ground, since in fact it is not possible to find a Noether 
symmetry and, therefore, an analytical exact general solution. This means that 
including radiation, so extending the model backwards, at least needs numerical 
techniques, if one also wants to retain the other two components, dust and scalar 
field. 
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Figure 1 : The value of u) must be <C 1 in order to obtain physical values for the 
scalar field density. 
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Figure 3: Evolution of the acceleration with time. The present time is to = 1. 
In the future the acceleration shows damped oscillation, with zero limit. 
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Figure 4: Scalar field density versus scale factor. The (black) thick dashed line 
is the matter density. The (blue) thin dashed line is drawn for reference and 
shows that the asymptotic value of the state equation for the scalar field is the 
same as matter. 
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Figure 5: The state equation of the scalar field versus the scale factor. Ap- 
parently, the present time (marked by the Y axis) is in the middle of a phase 
transition between w = — 1 and oscillations. The next figure clarifies the situa- 
tion. 
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Figure 6: Same as Fig. 5, but now several "present" values of f2o<£ are used, in 
the range 0.01 -^0.93. This shows that the plots are very similar and that the 
coincidence is only apparent. 
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